Gaps between jets in hadronic collisions 



O. Kepka/'2'E1 C. Marquet,3'[t| and C. Royon^^'S 



^ IPNP, Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic 
^Center for Particle Physics, Institute of Physics, Academy of Science, Prague, Czech Republic 
"^Physics Department, Theory Unit, CERN, 1211 Geneve 23, Switzerland 
I RFU/ Service de physique des particules, CEA/Saclay, 91191 Gif-sur-Yvette cedex, France 

We propose a model to describe difTractive events in liadron-hadron collisions where a rapidity gap 
is surrounded by two jets. The hard color-singlet object exchanged in the t-channel and responsible 
for the rapidity gap is described by the pQCD Balitsky-Fadin-Kuraev-Lipatov Pomeron, including 
corrections due to next-to-leading logarithms. We allow the rapidity gap to be smaller than the 
inter-jet rapidity interval, and the corresponding soft radiation is modeled using the HERWIG 
Monte Carlo. Our model is able to reproduce all Tevatron data, and allows to estimate the jet-gap- 
jet cross section at the LHC. 



A large effort has been devoted to understand the QCD dynamics of rapidity gaps in jet events since such processes 
were observed in p-l-p collisions at the Tevatron more than 10 years ago [l|, While describing diffractive processes 
in QCD has been a challenge for many years, the presence of a hard scale in so-called jet-gap-jet events for instance, 
brings hope that one could be able to understand these with perturbative methods. However, after many theoretical 
investigations, there is still no consensus on what the relevant QCD mechanism really is. 

In a hadron-hadron collision, a jet-gap-jet event features a large rapidity gap with a high-pT jet on each side 
(pt3>Aqc'd)- Across the gap, the object exchanged in the channel is color singlet and carries a large momentum 
transfer, and when the rapidity gap is sufficiently large the natural candidate in perturbative QCD is the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) Pomeron [3]. Of course the collision energy ^/s should be big {^/s ^ Et) in order 
for jets to be produced along with a large rapidity gap. Such events are expected to be produced copiously in p -\- p 
collisions at the LHC. 

To compute the jet-gap-jet process in the BFKL framework, one has first to address the problem of coupling the 
BFKL Pomeron to partons, as opposed to colorless particles. Indeed, BFKL calculations usually use the fact that 
impact factors, which describe the coupling of incoming and outgoing particles to the BFKL Pomeron, vanish when 
attached to gluons with no transverse momentum. This is a property of colorless impact factors. For instance, this 
is what allows to turn the Feynman-diagram calculation of the BFKL Pomeron into a conformal-invariant Green 
function Consequently, this BFKL Green function cannot be hooked to colored particles, and should be modified 
accordingly first. The Mueller- Tang (MT) prescription Q is widely used in the literature to couple the BFKL Pomeron 
to quarks and gluons. 

On the phenomenological side, the original parton-level MT calculation was not sufficient to describe the Tevatron 
data. A first attempt to improve it was proposed in , parton showering and hadronization were taken into account 
using the HERWIG Monte Carlo program 0. An agreement with data could only be obtained if the leading- 
logarithmic (LL) BFKL calculation was done with a fixed value of the coupling constant as, which is not satisfactory, 
as next-to-leading logarithmic (NLL) BFKL corrections are known to be important. In addition, only the leading 
conformal spin {p = 0) was taken into account. In [s!], it was shown that a good description of the data could be 
obtained when some NLL corrections were numerically taken into account in an effective way 0, but the full NLL- 
BFKL kernel [l^ could still not be implemented. As a result, these tests on the relevance of the BFKL dynamics 
were not conclusive. 

In the most recent phenomenological work on the subject [lH , the full NLL-BFKL kernel was implemented including 
all conformal spins, al ong with the collinear improvements necessary to remove spurious singularities and obtain 
meaningful results [l^, Il3| . However, the results of [ll| remained at parton level, therefore the fact that the rapidity 
interval between the jets can be larger than the rapidity gap could not be implemented. The DO measurement could 
nevertheless be reasonably well described, while the CDF data was not considered. The purpose of this letter is 
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FIG. 1: Production of two jets surrounding a large rapidity gap in a hadron-hadron collision, ^/s denotes the collision energy, 
PTi iVi) ^-nd PT2 iVi) the transverse momenta (rapidities) of the jets and xi and X2 are their longitudinal momentum fraction 
with respect to the incident hadrons. The rapidity interval between the jets Arjj is bigger than the rapidity gap A?7g. 



to improve the model by taking into account parton showering, hadronization effects and jet reconstruction, which 
is necessary to make more precise comparisons with data. We will interface the collinearly- improved NLL-BFKL 
parton-level results of [ll| with the HERWIG Monte Carlo, as was done in with the LL-BFKL calculation. We 
shall compare our results to both DO and CDF data. 

The plan of the letter is as follows. In section II, we recall the phenomenological NLL-BFKL formulation of the 
jet-gap-jet cross section and in section III, we explain how it is embedded into the HERWIG Monte Carlo program. In 
section IV, we present successful comparisons with all Tevatron data, which allow to fix the absolute normalization in 
our model. Predictions for the jet-gap-jet cross section at the LHC are presented in Section V. Section VI is devoted 
to conclusions and outlook. 



II. THE JET-GAP-JET CROSS SECTION IN THE BFKL FRAMEWORK 



The production of a rapidity gap between two outgoing jets in a hadron-hadron collision is pictured in Fig. [TJ with 
the different kinematic variables. We denote the collision energy, px^ and the transverse momenta of the 
two jets and xi and X2 their longitudinal fraction of momentum with respect to the incident hadrons. The rapidity 
interval between the two jets is Arjj — \n{xiX2s/pTiPT2)- At the parton level (see Fig.l in [IH), PTi = —PT2 — Pt, 
and the rapidity gap coincides with the rapidity interval Arj = la{xiX2S /p'^) between the outgoing partons that will 
initiate the jets. The hadronization of the partons into jets reduces the size of the rapidity gap to Arjg. 

In this section, we deal with the parton-level cross section 

dx,dx2dp\ = SUsAx,,VT)J.!f{^2,VT)^^. (1) 

where the functions fef f{x,p'^) are effective parton distributions that resum the leading logarithms log{p^ / Aq(-; jj) . 
They have the form 



feff{x,^l^) ^ g{x,^?) + J;^{q{x,^l'^) + q{x,^i^)) , (2) 



where g (respectively g, q) is the gluon (respectively quark, antiquark) distribution function in the incoming hadrons, 
and evolves according to DGLAP evolution [T^. Even though the process we consider involves moderate values of 
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xi and X2 and the perturbative scale ^ ^qcD' 'which we have chosen as the factorization scale, the cross section 
(H)) does not obey collinear factorization. This is due to possible secondary soft interactions between the colliding 
hadrons which can fill the rapidity gap. Therefore, in ([T|), the collinear factorization of the parton distributions /e// is 
corrected with the so-called gap-survival probability 5, which we assume depends only on ^/s as in standard diffractive 
calculations. Since the soft interactions happen on much longer time scales, the factor S is factorized from the hard 
part da^^^^^ /dp^. This hard cross section is given by 

'J^.J-lMA,.pif (3) 

in terms of the gg gg scattering amplitude A{Ar],p^). The two measured jets are initiated by the final-state gluons 
(or quarks), parton showering and hadronization effects will be discussed in the next section. 

In the following, we consider the high-energy limit in which the rapidity gap A77 is assumed to be very large. The 
BFKL framework allows to compute the gg — gg amplitude in this regime, and the result is known up to NLL 
accuracy. We note that there exist other QCD-based approaches to compute the jet-gap-jet cross section [15|. Let 
us first point out that in general collinear and fcr-factorization are two distinct schemes to factorize a hard process 
from a soft process (as is the case for the proton structure function F2), and should not be mixed. But the process 
we are investigating is different: collinear factorization is used to separate the hard part from the soft part, and 
/cT-factorization is only used within the hard part itself. It allows to factorize the amplitude A{Arj,p^) into three 
hard pieces: two impact factors defined order-by-order with respect to as, and the BFKL Green function where a 
resummation of leading (and next-leading) logarithms is performed. 

Since in our calculation the BFKL Pomeron is coupled to quarks or gluons, the BFKL Green function cannot 
be used as it is and should be modified. The transformation proposed in [5] is based on the fact that one should 
recover the analiticity of the Feynman diagrams. It was later argued that this prescripion corresponds to a deformed 
representation of the BFKL kernel that indeed could be coupled to colored particles and for which the bootstrap 
relation is fullfiled ^1^. Applying the MT prescription at NLL leads to 

2 167V^7ra|(p|) f dj b^-(7-l/2)^]exp{a(p|)xe//[2p,7,«(PT)]Ar?} 

^ Cfp\ [(7 -1/2)2 -(p- 1/2)2] [(7 -1/2)2 -(p+ 1/2)2] U 

with the complex integral run ning along the imaginary axis from 1/2— ioo to l/2-|-ioo, and with only even conformal 
spins contributing to the sum |17| . The running coupling is given by 



aip'r) = = [b\og{pl/KlcD)] , ^^N ' 



It is important to note that in formula we used the leading-order non- forward quark and gluon impact factors. 
We point out that the next-to-leading-order impact factors are known [ISJ , and that in principle a full NLL analysis 
is feasible, but this goes beyond the scope of our study. 

The NLL-BFKL effects are phenomenologically taken into account by the effective kernels Xeff{p, 7, ct)- For p = 0, 
the scheme-dependent NLL-BFKL kernels provided by the regularisation procedure XNLLij,!-^) depend on w, the 
Mellin variable conjugate to exp(A77). In each case, the NLL kernels obey a consistency condition [12l | which allows 
to reformulate the problem in terms of Xe//(7, a) (see also [H, [l^ for different approaches). The effective kernel 

M f (7 , a) is obtained from the NLL kernel XAfLL(7,^) by solving the implicit equation Xeff = Xnll{i,ol Xeff)- In 
|22| . the regularisation procedure has been extended to non-zero conformal spins and the kernel xmllIp, 1, ^) was 
obtained from the results of [2Q,]. The formulae needed to compute it can be found in the appendix of [2l| (in the 
present study we shall use the S4 scheme in which xnll is supplemented by an explicit a dependence, the results in 
the case of the S3 scheme are similar). Then the effective kernels Xe//(?5, 7, o;) are obtained from the NLL kernel by 
solving the implicit equation: 

Xeff = Xnll{p, 1, a Xeff) ■ (6) 

Similar NLL-BFKL phenomenological studies have been carried out with Mueller-Navelet jets in hadron-hadron 
collisions [2l|,[22], forward jet production in deep inelastic scattering ^23, 24], and the proton structure function 25,]. 
While in the F2 analysis the NLL corrections did not really improve the BFKL description, it was definitively the 
case in the forward-jet study. In the Mueller-Navelet jet case, NLL corrections dramatically change the predictions, 
even more so in the full calculation when NLO impact factors are also implemented ]^26]. In fact, these results cast 
strong doubts on the fact that Mueller-Navelet jets are a good observable to unambiguously observe BFKL effects, 
leaving the jet-gap-jet measurement as perhaps the new candidate. 
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In the LL-BFKL case that we consider for comparisons, the formula for the jet-gap-jet cross section is formahy the 
same as the NLL one, with the foUowing substitutions in 

Xeffip, 1, a) xll{p, 7) = 2-0(1) - -0 ^1 - 7 + - "0 (^7 + , "(fc^) ^ a = const, parameter, (7) 

where "0(7) = dlogT{'y)/d'y is the logarithmic derivative of the Gamma function. In this case, the coupling a is a 
priori a parameter. We choose to fix it to the value 0.16 obtained in [2^ by fitting the forward jet data from HERA. 
This unphysically small value of the coupling is indicative of the slower Bjorken-x dependence of the forward-jet data 
compared to the LL-BFKL cross section, when used with a reasonable a value. And in fact, the value a — 0.16 
mimics the slower energy dependence of NLL-BFKL cross section (in this case the average value of a is about 0.25), 
which in the forward-jet case is consistent with data. Therefore in both the LL- and NLL-BFKL cases, one deals with 
one-parameter formulae: the absolute normalization which is not under control. In the NLL case, this is due to the 
fact the we do not use NLO impact factors. 

Finally, to compute the cross section ([T]), we use CTEQ parton distribution functions [l^l, and we take S" = 0.1 for 
the gap-survival probability at the Tevatron and S = 0.03 at the LHC. More details on the parton-level computations 
can be found in fl^, such as the importance of the different conformal spins in (|4]), or the uncertainty due to the 
choice of the renormalization scale. In this work, the goal is to obtain hadron-level results by interfacing ([1]) with the 
HERWIG event generator. 



III. IMPLEMENTATION OF THE NLL-BFKL FORMULA IN HERWIG 



The parton-level calculation presented in the previous section leads by definition to a gap size Arj equal to the 
interval in rapidity between the partons that initiate the jets. At particle level, it is no longer true. Due to QCD 
radiation and hadronisation, the jets have a finite size, and the gap size Arjg is smaller than the difference in rapidity 
between the two jets Ar/j (see Fig. [Ij. This has an important consequence: to be able to compare the NLL-BFKL 
jet-gap-jet cross sections with CDF and DO measurements at the Tevatron, it is needed to embed our formulae in a 
Monte Carlo code. For instance, the DO collaboration selects events with a gap devoid of any activity in the [—1,1] 
region in rapidity while they require the jets to be separated by at least 4 units in rapidity. To take into account these 
effects, we implemented the NLL-BFKL cross section in the HERWIG Monte Carlo, and made our analysis at the 
particle level, after hadronisation. Since this procedure of going from parton-level to hadron-level is quite sensitive to 
the way jets are reconstructed, we use the same jet algorithm as experimentally used. 

Practically, in order to implement our formalism in HERWIG, we modified the HWHSNM function which imple- 
ments the matrix element squared for color-singlet parton-parton scattering 7]. Formula (jS]), which gives the BFKL 
da/dp^ cross section is too complicated to be implemented directly in HERWIG since it involves an integration in 
the complex plane over 7, and it would take too much computing time to generate many events. To avoid this issue, 
we parameterized da/dp^ as a function of the parton px and Ar] between both partons at generator level. Denoting 
^(Pt) ~ "^(Pt)^'?/^, the parametrization used is 

^ = "fjy^ [a + ^PT + + (d + epT + / VP?) X z + (g + hvT) x + (i -f 3^/P^) x 2^ + exp(fc lz)\ . (8) 

This formula is purely phenomenological, not motivated by theory, and was just introduced to obtain a very good 
while fitting ([8]) to the full expression of da/dp\. To perform the fit, 2330 points were used for parton pT ranging 
from 10 to 120 GeV, and A77 up to 10. The values of the different parameters were implemented in the HERWIG 
Monte Carlo. To summarize, we input into HERWIG the NLL-BFKL parton-level cross section which depends on pT 
and Ar/, and the output depends on Ar^g, Ar/j, and the jets transverse momenta pt^ and pji- Further integrations of 
these kinematic variables are performed to obtain the different observables discussed in the next section, taking into 
account experimental cuts. 



IV. COMPARAISON WITH TEVATRON DATA 



The DO collaboration has performed a measurement of the jet-gap-jet event ratio, defined as the ratio of the jet- 
gap-jet cross section to the inclusive di-jet cross section, as a function of the transverse energy of the second-leading 
jet, that we denote Et, and also as a function of the rapidity difference Arjj between the two leading jets At 
least two jets are reconstructed in the DO calorimeter with Et > 15 GeV for the second leading jet. In addition, 
the two jets are required to be in the forward regions and in opposite hemispheres, by requesting 1.9 < |?/i,2| < 4.1 
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FIG. 2: Comparisons between the DO (left) and CDF (right) measurements of the jet-gap-jet event ratio with the NLL- and 
LL-BFKL calculations. The NLL calculation is in fair agreement with the data while the LL one leads to a worse description. 



and r]iri2 < 0. The difference in rapidity between both jets Aryj was imposed to be larger than 4, and a rapidity gap 
between at least rj ~ —1 and 77 = 1 is required. The data are presented as a function of the second- leading-jet Et, or 
as a function of Aj]j in which case a Iow-Et and a high- jet samples were used (low Et means 15 < Et < 25 GeV 
and high Et means Et > 30 GeV, those cuts applying to both jets). 

To compare directly the DO measurement with the NLL-BFKL calculation implemented in HERWIG, we compute 
the following ratio 

_ NLL BFKL Herwig ^ LO QCD 

~ Jet Herwig ^ NLO QCD ' ^' 

where NLL BFKL Herwig and Dijet Herwig are the jet-gap-jet and the inclusive di-jet cross sections obtained with 
HERWIG, respectively. To take into account NLO QCD effects, we also correct the ratio R by the LO/NLO QCD 
di-jet cross section ratio obtained with the NLOJet-f+ program (2^. The same method applies for the LL-BFKL 
cross section calculation. The comparison between our calculations and the DO data is given in the left plot of Fig. [21 
after the overall normalization was adjusted in the two cases (LL and NLL). We find that there is a good agreement 
between the NLL-BFKL calculation and the data whereas the LL-BFKL calculation leads to an Et dependence which 
is too flat. Moreover, in the LL case it is difficult to accomodate all the data with a single overall normalization factor, 
while in the NLL case this is not a problem. 

The comparison with CDF data is given in the right plot of Figl^l in this case the normalization of the different 
data sets is arbitrary. The CDF collaboration also measured the jet-gap-jet cross section requesting a gap between -1 
and 1 in rapidity, but used a higher jet Et threshold of 20 GeV and a lower acceptance in jet rapidity between 1.8 to 
3.5, compared to the DO measurement. The CDF requirement on the minimum rapidity interval between the two jets 
is then Arjj > 3.6, compared to 4 in the case of DO. CDF measured the jet-gap-jet event ratio as a function of the 
average jet transverse momentum (this is what Et denotes in the case of CDF), the average transverse momentum 

of the third jet Ej, when there is one in the event, and also as a function of Arjj. The conclusion remains the same 
as in the case of DO data, namely that the NLL-BFKL formalism leads to a better description than the LL one. 
However, it is worth noticing that we are not able to describe the full Arjj dependence and especially the decrease 
at high A?7j, which is somehow in disagreement with the DO measurement. Further measurements in progress in the 
CDF collaboration will be very useful to understand these differences. Parton showering and hadronization effects 
are crucial in order to obtain this level of agreement with data. 
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FIG. 3: Predictions of our model for the ratio of the jet-gap-jet to the inclusive-jet cross section at the LHC, as a function of 
the second-leading-jet transverse energy Et (left), and of the rapidity difference between the two leading jets Arjj (right). 



V. PREDICTIONS FOR THE LHC 



Using the normalizations obtained from the fits to the DO data, we are able to predict the jet-gap-jet event ratio at 
the LHC. In doing so, we also take into account the fact that the gap-survival probability is smaller by a factor 10/3 
(this is estimated for a collision energy of y/s = 14 TeV). Requesting both jets to have Et > 20 GeV, and the jet 
rapidities to obey 2 < \r]i^2\ < 5 and r]ir]2 < 0, the values of jet-gap-jet event ratios are shown in Fig. [3l as a function 
of Et for different Ai]j ranges (left plot), and as a function of Arjj for different Et ranges (right plot). The main 
feature of the predictions of our model is that jet-gap-jet event ratio is about 0.002 and does not vary a lot with Et 
or A77,/, except for the LL-BFKL predictions which increase a bit with Aryj. By contrast, the parton- level predictions 
obtained in [ll| featured an increase of the ratio with both Et or Ar]j. 

Let us add a word of caution about these predictions. We have assumed a specific value for the gap survival 
probability at the LHC (namely 0.03), however this value shows large theoretical uncertainties [1^. It will be 
measured by the LHC experiments and our cross section predictions will have to be modified once such measurements 
are performed. For instance, our cross section should be reduced by a factor 10 if the survival probability is found to 
be ten times smaller than 0.03 for a center-of-mass energy of 14 TeV. In addition, our prediction assumes that the 
inclusive jet cross section at the LHC can be correctly described by the Herwig Monte Carlo, as is the case at the 
Tevatron [s^. This can be again tested at the LHC using the first data. One first indication for a center-of-mass 
energy of 7 TeV was given by the measurement of the difference in azimuthal angle in dijet events [3l| which seems 
to be well described by the Herwig Monte Carlo. However, the same measurement as well as cross section comparison 
between data and MC need to be done at a higher center-of-mass energy at 14 TeV to make sure that the Herwig 
Monte Carlo can describe the inclusive jet measurement. If some discrepancy need to be accounted for, this has to 
be taken into account as well in our prediction of the jet gap jet cross section ratio. 



VI. CONCLUSIONS 

We have embedded the parton-level NLL-BFKL calculation of [llj into the HERWIG Monte Carlo program, in 
order to obtain hadron-level results for the jet-gap-jet cross-section in hadron-hadron collisions, corresponding to 
the production of two high-p^ jets around a large rapidity gap. The NLL-BFKL effects are implemented through a 
renormalization-group improved kernel in the S4 scheme, while the Mueller- Tang prescription is used to couple the 
BFKL Pomeron to colored partons, described with only LO impact factors. 
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After adjusting one parameter, the overall normalization, the NLL-BFKL calculation is able to describe all Tevatron 
data, except the higher end of the Ajyj dependence measured by CDF. This still provides an improvement compared 
to the LL-BFKL calculation (obtained with the fixed value of the coupling a = 0.16), which in addition features an 
Et dependence that is too flat compared to the DO data. We presented predictions which could be tested at the LHC, 
for the same jet-gap-jet event ratio measured at the Tevatron, but for larger rapidity gaps. 

We noticed that going from parton-level to hadron-level is necessary in order to obtain a global description of the 
Tevatron data, hence our results supersede those of Ref. fTT'|. Considering LHC predictions, hadron-level calculations 
show almost no dependence of the jet-gap-jet event ratio (which is about 0.002 in the kinematic range we considered, 
and for the experimental cuts we used) , while parton-level calculations showed an increase of this ratio with both Et 
and Arjj. This should provide a strong test of the BFKL regime. 
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